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^ ■ Abstract: SPM (Sand Pile Model) is a simple discrete dynamical system used in 

physics to represent granular objects. It is deeply related to integer partitions, and 
many other combinatorics problems, such as tilings or rewriting systems. The evolution 
of the system started with n stacked grains generates a lattice, denoted by SPM{n). 
We study here the structure of this lattice. We first explain how it can be constructed, 
by showing its strong self-similarity property. Then, we define SPM{oo), a natural 
extension of SPM when one starts with an infinite number of grains. Again, we give 
an efficient construction algorithm and a coding of this lattice using a self-similar tree. 
. The two approaches give different recursive formulae for \SPM{n)\. 

Keyvirords: SPM, Sand Pile Model, Lattice, Integers partitions, CFG, Discrete 
Dynamical Systems. 
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O 1 Introduction 

O , 

95 ! 1.1 Motivations and context 

O 

■ In 1987, Bak, Tang and Wiesenfeld [BTW87] introduced the important notion of self- 

\ organisation criticality (soc): when certain systems in a steady state (named critical 

C/3 ■ state) are slightly perturbated, they evolve back to another steady state. This evolution 

implies some arbitrarily high modifications of the system. 

The typical example is an avalanche on a sand pile. At first, the pile is in a steady 
k> \ state and the perturbation consists in adding a grain on the pile. As a consequence, 

j_j ■ the pile evolves to a new steady state, with an avalanche starting where the grain 

d I was dropped. The fact that this avalanche size may be arbitrarily high is the main 

characteristic of SOC systems. 

Since the appearance of this paper, many physicists and biologists have recognized 
these properties in natural systems, and the SOC family still grows (see [Jen98] , [Tan93] 
for example) ; many publications on this topic appeared recently [Jen98] , [Dur97] , [Bak97] , 
[Tur97]. These phenomena are of particular interest in surface grow studies [BS95], in 
geophysics [BT98], in plasma confinment, in astrophysics, and many other, including, of 
course, studies of granular systems like dunes [Dur97] and molecule agregation [BS95] . 

The essence of these phenomena is captured by a well known model in game theory 
and combinatorics, the Chip Firing Game (cfg). The most general notion of CFG 
is a directed graph G = {V,E) where a threshold Sy and a load l{v) are given to each 
vertex v. Intuitively, l{v) represents the number of chips stored at v. The game evolves 
with respect to the following rule: if v Cz V contains more than Sy chips, then it gives 
of them to its neighbours, i.e. the load of the vertex v is decreased by and the load 
of each of its neighbours is increased by ^ where ny is the number of v neighbours. In 
general, one takes Sy — ny, but 6y ^ oo if ny — {v is then called a sink). See Figure 1 
for an example of such a CFG. 
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Figure 1: An example of CFG. Each vertex load is written in this vertex. 



Under certains conditions, the CFG converges to a steady state (see for example 

[Eri93]). The addition of one chip on a vertex w 6 V when the system is in a steady state 
causes a redistribution of the chips. During this redistribution, an arbitrary number of 
vertexs may be concerned. To give an example, we can consider the case where /(z/) =0 
and l{v) = 5v y V £ V \ {v}. If one adds successively 5^, ~ 1 chips on i/, the only 
concerned vertex is v (the system remains steady). If one adds one more grain, every 
vertex will be concerned (if G is connected). Such a diffusion can be arbitrary large 
[GM97] , depending on the initial state of the system, and is always started by addition 
of one grain. Such a propagation is called an avalanche. 

A particular case of this model is widely studied: the sand pile model on a rectangular 
grid^. The graph G in this case is undirected. It is a rectangular finite lattice and the 
value of 6v and of riy is 4 for all vertex v except one singular vertex v which is linked 
once to any vertex on the border of the lattice and twice to the four corners, and such 
that dy = oo. The distinguished vertex acts like a sink: it never gives away any of 
its grains and could be considered as collecting the grains that leaves the system. If 
the load of a vertex inside the lattice is more than 4, then it gives one grain to each 
of its four neighbours (sec Figure 2 for an example where the distinguished vertex is 
not represented since its load does not influence the evolution of the system). This 
is the model deeply studied by Dhar [DM90], [DRSV95]. In particular, one can show 
that adding a grain turns the system into an unsteady state, and that after auto- 
reorganisation it reaches a new steady state [Mar92] . This confirms that we are in the 
SOC context. Cori and Rossin [CR98] generalized this notion to any rooted graph and 
obtained similar results. 
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Figure 2: Example of (parallel) sand pile. 

Another special case of cfg is the Sand Piles Model (spm). The graph G in this 
case is an undirected chain, infinite on the right: V = N, E = + 1) V « G N}, 
Sy = 2 for all w > and Sq = 00 (see Figure 3(a)). This model is equivalent to the 
following. Consider an infinite chain of columns, each containing a vertical pile of 
grains. The height difference between the column Cj and its right neighbour column 
Cj+i is denoted by d{i). If d{i) is greater than or equal to 2 then a grain falls down 
from d to Cj+i (see Figure 3(b)). If i > 0, we call Cj-i the left neighbour of c, and then 
d{ci-i) and d{ci+i) are increased by 1 while d{ci) is decreased by 2. We find again our 
initial definition of SPM, with a coding of the pile by height differences. 

Notice that spm is less general than the lattice sand pile but is not a particular case 
of it: the number of vertices v with l{v) 7^ is not bounded in SPM, and it effectively 
grows with the number of grains. Moreover, the model SPM has no sink, which is 

^The standard term is lattice but, since we will use orders theory in the following, where the word 
lattice takes another meaning, we use here the word grid. 
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a fundamental difference. If wc represent a pile by the t-uple of its columns height, 
each configuration of the pile represents a partition of the total number of grains (see 
Figure 3(c)). 

(a) 0| 2| 1| 1| OTOP ---^ 1| 0| 2| 1| 0]l)P ---- ni| 1| 0| 2| 0| OP ---- n'F'l 'I "I 'I "I 



(b) 




W (4,4,2,1) (4.3.3.1) » (4.3.2.2) » (4,3,2,1,1) 



Figure 3: Three ways to see spm. 

In computer science, the cfg models several problems and is applied in several algo- 
rithms (see for example [GP97]). spm itself admits natural interpretations in algorithmic 

terms. Wc give here two examples about dynamical distribution of jobs on a proces- 
sors network [Hua93] [DKTR95] [GMP98a]. Each column of a sand pile represents a 
processor, a grain represents a job. One can imagine the processors are connected on 
a ring (like Token Ring): each processor can only communicate directly with its right 
neighbour. It corresponds to the move of a grain from one column to another. Since 
only neighbour processors can communicate, the communications can be processed in 
parallel and the parallel SPM is a good model for this problem [DL96]. If on the contrary 
the communication medium is a shared bus (like Ethernet or certain multiprocessors), 
we can study the evolution of sequential spm to avoid collisions. 

In the following, we are going to discuss some lattice properties of the above dy- 
namical systems. Let us recall that a lattice can be described as a partial order such 
that two elements a and b admit a least upper bound (called supremum of a and b 
and denoted by sup(a, b)) and a greatest lower bound (called infimum of a and b and 
denoted by inf (a, b)). The element sup(a, 6) is the smallest element among the elements 
greater than both a and b. The element inf (a, 6) is defined similarly. A useful result 
about finite lattices is that a partial order is a lattice if and only if it admits a greatest 
element, and any two elements admit a greatest lower bound. For more details, see 
for example [DP90] . The fact that a dynamical system has the lattice property implies 
some important properties, such as convergence. 



1.2 Our model: known results 

Our model is the standard sequential SPM; it consists of an infinite number of ordered 
columns, each containing a certain number of grains. Only the first k columns are non- 
empty, so the state of the system is described by the fc-uple s = (si, S2, . . . , Sfe) where 
Si is the number of grains in the column i for 1 < i < k. 

The system is initially in the state N = (n). This means that all the grains are in 
the first column. At each step, the system evolves with respect to the following rule: 
one grain can fall down from column i to column i -|- 1 if and only if Sj+i — Sj > 2. This 
rule defines a covering relation on the set of reachable configurations. The refiexive and 
transitive closure of this relation is an order, called the dominance order [GMP98b]. 
The set of reachable configurations from the partition (n) with this order is then a 
lattice denoted by SPM{n) [GMP98b]. 

Let s = (si, . . . Sfe) be a sand pile, the height difference of s at i, denoted by di{s), 
is the integer Sj — Sj+i (with the assumption that Sk+i = 0). We will say that s has a 
step (resp. plateau, resp. cliff) at i if and only if its height difference at i is 1 (resp. 0, 
resp. > 2). We extend these definitions by saying that s has stairs (resp. a plateau) 
at the interval if and only if s has a step (resp. plateau) at k for all i < k < j. 
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The integer j — i + 1 is called the length of the stairs (resp. plateau). See Figure 4 for 
examples. 



stairs plateau I t I t I 
cliff plateau step plateau cliff 



Figure 4: Examples for the definitions. 

The evolution rule of a sand pile s = (si, . . . , Sj, Sj+i, . . . , Sk) is then: one grain can 
fall from one column to the column on its right if and only if it is at the top of a cliff. 

Such a transition is denoted by — ^ where i is the number of the column from which 
the grain falls. The sand pile s' is called a successor of s, and Succ{s) denotes the set 
of all successors of s: 

Succ{s) = {s'\s^s',ieN}. 
See Figure 5 for an example. 




Figure 5: Transitions, successors. 

Let us now introduce a few notations. If s is a partition of n then s^* is the partition 

of n + 1 obtained by adding a grain on the ?'-th column of s (if it is a partition). In 
other words, if s = (si, . . . , Sj, . . . , Sk) then s^' = (si, . . . , s, + 1, . . . , Sfe). We also define 
gli — g where S is a set of partitions. 

We will denote by e(s) the largest integer such that s has stairs at [1, e(s)]. We also 
define Pi as the set of the sand piles in SPM{n) that begins with stairs of length (at 
least) i. In other words, = {s g SPM{n) \ e{s) > i}. 

Charaterisations of the fixed point of the system, the minimum element of the cor- 
responding lattice, and of its elements arc also known: 

Theorem 1 ([GK93]) The fixed point of SPM{n) is: 

So = {k,k - 1, ... ,p+ l,p,p,p - 1, ... ,2,1) 

where k is the maximal integer such that Sq is a sand pile of n grains, i.e. k is the 

integer such that ^''^^'^^ < n < Lt+Dlt+^ _ 

Theorem 2 ([GMP98b]) A partition s belongs to SPM{n) if and only if: 

• s does not contain any sequence p,p,p or p,p, {p — 1), {p — 1). 

• there is at least one cliff between two consecutive sequences p,p and q,q. 
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In this paper, we study the structure of SPM{n). In particular, we show in the 
next section how SPM{n + 1) can be constructed from SPM{n), thus we obtain an 
algorithm that constructs SPM{n) for any integer n. Afterwards, wc define; a natural 
infinite extension, SPM (oo), when the system is started with an infinite column of 
grains. The study of the structure of SPM{oo) permits more remarks on the self- 
similarity of the set. During this study, we obtain interesting recursive formula for 
\SPM{n)\. 



2 From SPM{n) to SPM{n + 1) 

The goal of this section is the construction of the lattice SPM{n + 1) from SPM{n). 
We will construct the graph of the transitive reduction of the lattice, i.e. the graph of 
its order relation, without the reflexive edges {x — > x) and the transitive ones {x — >■ z 
when X — > y and y — > z). Each edge of this graph is equivalent to a transition of 

the SPM system. Therefore, we will label the edge s — ^ s' with the number i of the 
column of s from which the grain falls in order to obtain s' . We will call the obtained 
labelled graph the diagram of the lattice. Wc first give some preliminary results, then 
we notice that SPM(n) is a good starting point to construct SPM{n + 1), and we give 
a method to obtain SPM{n + 1) from SPM{n). Finally, we inspect more deeply the 
construction algorithm and show a strong self-similarity in each lattice SPM(n). This 
similarity induces a first recursive formula for the cardinality of SPM(n). 



2.1 Preliminaries 



Let us study what happens when wc add one grain on the i-th column of a sand pile 
s = (si, . . . , Sj, . . . , Sfe) such that e(s) >i — l. We obtain the sand pile s-l-* = (si, . . . , Si + 
1, . . . , Sfc). We want to determine all the possible transitions from this partition, knowing 
the possible ones from s. Three cases are possible (as shown in Figure 6) corresponding 
to the three following propositions. Recall that we only consider sand piles s with 
e(s) > i — 1, since it will be the case of interest for the rest of the paper. 
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Figure 6: The three cases considered. 



Proposition 1 (plateau) Let s e SPM{n) such that e(s) >i — l. If s has a plateau 

at i then the possible transitions from s-'-' are the same as the possible transitions from 

s. Moreover, if s t then t^K In other words, Succ{s^*) = {Succ{s))^' and 

the corresponding edges of the diagrams have the same labels. 

Proof: A transition is only possible if there is a cliff at the column i. Now, the 
set of the columns where s has a cliff is equal to the set of the columns where s ^^ has a 
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Proposition 2 (cliff) Let s G SPM{n) such that e(s) > i — 1. If s has a cliff at i 
then 

1. The possible transitions from s^' are the same as the possible transitions from s 

and if s t then s^* t^K In other words, Succ{s^*) = {Succ{s))^* and the 
corresponding edges of the diagrams have the same labels. 

2. Moreover, if s^' s' then s^^ ~ — ^ — ^ ' ^ s' . In other words, s' is reachable 
in SPM(nY^ from s^'^ via a path labelled i ■ i — 1 ■ ... •2-1. 

Proof: 

1. A transition — ^ is only possible if there is a cliff at the column i. Now, the set 
of the columns where s has a cliff is equal to the set of the columns where s^^ has 
a cliff. 

2. The sand pile s' is equal to s^^+^ and by hypothesis s has a cliff at i and stairs at 
[1, i — 1]. Therefore, in the sand pile s^^ a grain can fall from the column i. With 

this grain's fall, we create a cliff at i — 1, therefore a new grain can now fall from 
column i — 2. This process can be iterated to obtain s' at the end. To sum up, 
we can write: 

S^^ = (.Sl + 1, . . . , Sj-l, Sj,Si+l, . . . , Sfe) 

(si + 1, . . . , Si-i,Si - 1, Si+i + 1, . . . , Sfe) 

i— 1 2 

> . . . > (Si + 1, S2 - 1, . . . , Si, Si+i + 1, . . . , Sfe) 

(Si + 1 - 1, S2, . . . , S^,S^+l + 1, . . . , Sfc) = S^*+^ = s' 

It is obvious that all the partitions on this path belong to SPM{n)^'^ . 

□ 

Proposition 3 (step) Let s G SPM{n) such that e{s) > i — If s has a step at i then 
the possible transitions from s^* are the same as from s with an additionnal transition 
on the column i: s^^ s-l-'+i. 

Proof: The set of columns where s^' has a cliff is equal to the union of {i} and the set 
of columns where s has a cliff. □ 

2.2 Construction 

Using the preliminary results from the previous section, we will here obtain an algorithm 

for the construction of SPM{n + 1) from SPM{n). Wc first show that SPM{n)^^ is a 
good starting point for the construction of SPM(n + 1). Recall that SPM{n)^'^ is the 
set of partitions obtained by addition of one grain on the first column of each partition in 
SPM{n). Afterwards, we will use the previous propositions to add the missing elements 
and transitions in order to complete SPM{nY^ into SPM{n+ 1). 

Proposition 4 SPM{n)^^ is a sublattice of SPM{n+l). 

Proof: Let us recall that if a and 6 are two partitions of SPM{n) for a given n, then 
inf (a, b) is their first common descendant and sup(a, b) is their first common ancestor. 
To prove the claim, we must show that: 

• If inf (a, 6) = c is in SPM{n) then mi{a^\b^^) = c^i is in SPM{n+l). 
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• If sup(a, 5) = c is in SPM{n) then sup(a^\ 6^1) = c^^ is in SPM{n + 1). 



Recall that [GMP98b]: 



inf (a, 6) = c iff for all j one has 

i=l 



i=l i=l 



This implies that: 




+ 1) for all j, 



i—1 i=l i=l 



i.e. c^i is in SPM{n+l). 

Let now c be equal to sup(a, 6) (in SPM{n)) and d be equal to sup(aJ^\6^^) (in 
SPM{nY^). Wc will show that d = c^^ . We have c > a and c > 6, therefore c^^ > a^^ 
and cl-i > h^^ , which implies that c^^ > d. Let us begin by showing that di = ci + 1. 
We can suppose a\ > b\. The partition {a\,a\, ai — 1, oi — 2, . . . ) is greater than a and 
h, hence it is greater than c. This implies that ci = ai. Since a'^^ < li < c^^, wc have 
di = ai + 1. Let e = (di — 1, d2, d^, . . . ). Since d < c^^ , e verifies the characterisation of 
Theorem 2. Moreover, d > a^^ and d>b^^, hence e> a and e > b. This implies that 
e > sup(a, b) = c and that d > c-I^S and so = c-l^i . □ 

It is straightforward that each clement s of SPM{n + 1) is reachable from an clement 
of SPM(nY^ . Indeed, s is at least reachable from (n)^^ = (n+ 1). This shows that one 
can start the construction of SPM{n + 1) with SPM{n)^^ and then add the missing 
elements (see Figure 7 for an example). 

The construction procedure starts with the lattice SPM{n)^^ given by its diagram. 
Then, we look for those elements in SPM{n)^^ that have a successor out of SPM{n)^^ . 
The set of these elements will be denoted by Ji, with Ii C SPM{nY^ . At this point, we 
add all the missing successors of the elements of Ji . The set of these new elements will 
be denoted by C\ . Now, we look for the elements in C\ that have a successor out of the 
constructed set. The set of these elements is denoted by 12- We add the new elements 
(their set is denoted by C2), and we iterate this process until the set Ii is empty. 

More explicitly, in the i-th step of the procedure we look for the elements in Cj_i 
with missing successors and call Ii the set of these elements. We add the new successors 
of the elements of Ii and call the set of these new elements Ci. At each step, when 
we add a new element, we also add its covering relations. SPM{n +1) is a finite 
set, therefore this procedure terminates. At the end, we have obtained the whole set 
SPM(n + 1) with its order relation. 

Now, let us show how this completion of SPM{n)^^ to obtain SPM{n + 1) is im- 
plemented. Recall that any element t of SPM{n)^'^ is obtained from an element s of 
SPM{n) by adding a new grain on the first column. Three cases are possible: 

• s begins with a plateau or a cliff. Then, according to Proposition 1, the possible 
transitions from s^^ are the same as the possible transitions from s, and the suc- 
cessors of s-^^ are obtained by an application of to the successors of s. Therefore 
they arc already in SPM{nY^ (sec Figure 8(a)(b)). 

• s begins with a step. In this case, s is in Pi = SPM{n)di=i- Then, according 
to Proposition 3, only one successor of s^'^ is not yet in SPM{nY^: the successor 

obtained by a transition — This element is s^^ (see Figure 8(c)). It follows that 



h=Pl^ andCi =Pl^^ 
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4,3,2,1,1 



Figure 7: SPM{10y^ in SPM{U). 

This means that the first step of the construction consists in adding all the elements 
of Ci = . Notice that this set is added with a duplication of the order structure of 
P/^ = Ii. Indeed, it is clear that: 

(s^i i^i) iff (s^2 t^^) for aU s, i in Pi, and for all j. 

The following step consists in adding the missing successors of the elements of Ci and 
the missing transitions originating from them. The analysis of the three cases (plateau, 
cliff, step) shows that the only elements of Ci that do not have all their successors and 
transitions are: 

• The elements s^^ G Ci such that s has a cliff at 2. Indeed, such a s^^ does have 
all its successors in the lattice, but one transition is missing: the one labelled with 
2. In this case Proposition 2 shows that s^^ t where t is also obtained by 
s^^ t. Therefore, we have to add an edge from s^^ to an element t which 
is already in the lattice. We will call back edge such an edge. 
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(a) (b) (c) 

Figure 8: Successors of s^^ if s begins with (a) a plateau (b) a cliff and (c) a step: 
there is a new element only in this last case. 

• The elements s^^ G Ci such that s has a step at 2 (i.e. s begins with stairs of length 
at least 2 and hence is in P2^). According to Proposition 3, only one successor 
of each of these elements is not yet in the lattice: the successor obtained by the 
transition on the second column, i.e. the clement .s^''. Therefore, to complete the 
second step, we have to add the set Pj^^ (with the same order structure as Pj^^ ) to 
the existing lattice and connect the lattice to this new part by all the transitions: 

gl2 foj. all _,i2 g pi-Z _ 

This means that I2 = and C2 = ■ In general, the i-th step consists in adding 
the missing successors of the elements added at step i — 1 and the missing transitions 
originating from them. We show that the observed behaviour for the second step is 
general, and so the sets and C, can be characterized. 

Theorem 3 For all integer i, we have Ii = P/* and Ci = ph+^ . 

Proof: By induction: 

• The case i = 1 has already been studied. Notice that every covering relation 
concerning the new elements is of the following form: s^'^ — ^ s^^ where s € Pi. 

• Suppose the result is true for i — 1. We show that it is true for i. Consider 
Ci-i = Pill- Using Propositions 1, 2 and 3, we look for the successors of s^', 
with s e Pi-i- Three cases are to be considered: s can have a plateau, a cliff or 
a step at i. 

Plateau According to Proposition 1, Succ{s'^^) C ph^ = d^i. So, .s has no new 

successor, and ,s ^ Moreover, s t if and only if s^^ f}'' . 

Cliff According to Proposition 2, Succ{s^') C p/_:^ U SPM{n)^K So, s ^ h. 
Moreover, the edges of the covering relation originating from s^* are the 

same than the ones originating from s plus an additional one: s^^ — ^ s^*+^ . 
The element s^*+'^ is in SPM{nY^ according to Proposition 2. 
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stair According to Proposition 3, s^^ has a new successor, s^'+^, hence s^' € /». 
Moreover, the edges of the covering relation originating from s^' are the 

same as the ones originating from s plus an additional one: s^-* — ^ s^'+^. 
Therefore the element gJ^'+i is in Cj. 

Prom these three cases, we deduce the claim. q 

Wc have obtained a characterization of the sets and C'i . It is now straightforward 
that Algorithm 1 constructs the lattice SPM{n + 1) from SPM{n). Notice that we 
can obtain SPM{n) for an arbitrary integer n by starting from SPM{0) and iterating 
this algorithm. In the next sections, we will give more details about this construction. 
We will show that the complexity of Algorithm 1 is linear with respect to the number 
of newly added elements, and hence we have an algorithm that constructs SPM{n) in 
linear time linear with respect to \SPM{n)\. 

Algorithm 1 Incremental construction 



Input: SPM{n) 
Output: SPM{n- 
begin 



1) 



Pi 



while 7 7^ do 



add C with its covering relation; 
for each s-^' in I do 
|_ add the edge: s^^ — ^ s^*+^ e C 

for each s^^+^ in C s.t. di+i{s) > 2 do 
|_ add the back edge s^'+i ^ s' € SPM{nY^ 



i + 1; 



end 



2.3 Structure of the parts 

We will now study more deeply the construction procedure given above. We will obtain 
results on the structure of the Pi parts, which play an important role, and a recursive 
formula for \SPM(n)\. However, the results presented here are not necessary to under- 
stand the infinite extension presented in the second part of the paper. Therefore, the 
rest of this section can be ignored if the reader is mostly interested in the second part 
of the paper. 

In the previous section, we characterized the sets and C,. More can be said about 
the structure of these sets. In fact, since /j = P/* and Cj = we only have to 

study the sets Pi. We will show that these sets are disjoint unions of lattices, and that 
each of these lattices is obtained from a generating partition by iteration of the SPM 
rule. We will give the explicit characterisation of these generating partitions, as well as 
their number. 

Proposition 5 Pi is a disjoint union of lattices. 
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Proof: Let denote the set of all the elements of SPM{n) whose first two parts 
are k and k — 1. This is a non-empty subset of Pi. It is clear that if k k' then 
Qi,k n Qi^k' = 0, so Pi is the disjoint miion of the sots Qi,fe. 

Since Pi C SPM{n), the elements of Qi.fc verify the characterisation of Theorem 2; 
this implies that the maximal element g of Q\^k has the form 

g = {k,k — l,k — l,k — 2, . . . ,k — l,r) 

with I maximal (i.e. I such that r < k — I — 1), and k + {k — 1) + {k — 1) + {k — 2) + 
. . . + {k — I) + r = n. Then g belongs to SPM{n) if and only if k satisfies: 

fc-|-fc-l<n<fc + (fc-l) + (fc-l) + (fc-2) + ...-|-l 

which is equivalent to2fc-l<n< 

Let k be such an integer. Let us study the structure of Qi^k by considering its 
maximal element g, described as above. Let s be an element of Qi,k- It is clear that 
the prefix sums of s are less than or equal to the ones of g, so, according to [GMP98b], 
s can be obtained from g by the spm rule. Therefore, Qi.fc is the set of the elements 
of SPM{n) which can be obtained from g and whose first two parts are k and fc — 1. 
In other words, Qi^k is the set of the partitions reached from g by paths without any 
transition labelled 1 or 2. The element g is called the generating partition of Qi.fc. 
Therefore Qi.k is isomorphic to the lattice of the partitions of n — 2A; — 1 obtained from 
((73, . . . , (/„) by iteration of the SPM rule, and, in particular, Qi^k is a lattice. q 

More generally, let us denote by Qi,k the set Qi,k H Pi. 

Proposition 6 The sets Qi^k o.fe lattices with all transitions labelled with integers 
greater than k. Moreover, for all i, Pi is the disjoint union of the lattices Qi^k- 

Proof: Recall that Pj is the subset of Pi containing the partitions that begin with 

stairs of length z, and that Qi^k is the subset of Pi containing the partitions that begin 
with k,k — 1. So, Qi^k is the subset of Pi containing the partitions that begin with the 
stairs k,k — 1, . . . ,k — i. Therefore the maximal element of Qi^k has the form: 

g = {k, k — 1, . . . , k — i, k — i, . . . ,k — l,r) 

where I > i and I is maximal (i.e. r < k — I — 1) (we use here the same argument as in 
the proof of Proposition 5). Every element s of Qi.fe is reachable from this element g, 
and only transitions with labels greater than i arc needed to obtain s from g. Therefore, 
Qi^k is a lattice isomorphic to the lattice of partitions obtained from (k — i,... ,k — l,r). 
We have 

k 

where U denotes the disjoint union, and the Qi^k are pairwise disjoint, then: 

Pi = p^nPi = [_\{Qi,k n Pi) = □ Qi,k, 

k k 

and obviously the sets Qi^k are also pairwise disjoint for a fixed i. q 

We sum up these results in Figure 9. An example is given for n = 10 in Figure 10. 

We have defined the generating partition g of a set Qi,fe as the maximal clement of 
Qi^k- Therefore Qi^k is the lattice of the sand piles reachable from g by iteration of the 
SPM rule on the coulmns at the right of the second column. Using the characterisation 
of SPM partitions, we can now enumerate all the generating partitions for a given n. 
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Figure 9: Structure of SPM(n+l). 

Proposition 7 (Generating partitions) The number of generating partitions in the 
set SPM{n) m [f + 2 - ^J^2n\ . 

Proof: As seen above, the generating partitions in SPM(n) have the form: {k, k — 
l,k — l,k — 2, . . . ,k — l,r) for some I > and r < k — I — 1. Such integers k must verify 
k+{k — l)<n and k — l + '^^'^^^^ > n. Moreover, any of these k effectively corresponds 
to a generating partition. Therefore we have as many generating partitions as solutions 
to the system: 

\ k + {k-l)<n, 

thatis-| + yfT2^<A;< □ 

Wo have already seen that each Qij is a lattice and contains the lattices Qij for 
i > 1. The lattices Qij also verify: Qi+ij C Q^ j for all j. We will show how the Qij 
are generated during the construction. See Figure 11. 

In order to study the parts Pi and Qi j when n varies, let us extend our notations. 
We will denote by Pi{n) the parts Pi of SPM{n) {Pi{n) is the set of all sand piles with 
n grains that begins with stairs of length at least i). Likewise, Qij{n) denotes the part 
Qij of SPM{n). We can make three remarks: 

• The elements obtained from P,(n) by applying the operator any number of 
times do not belong to a Pi{m) with m > n since they begin with a cliff. 
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Figure 10: Example with n=10. 



• The elements of {Pi{n))^^ begin with a plateau at column 1 followed by a cliff 
at column 2. Then, if we apply to these elements we obtain sand piles which 
begin with stairs of length exactly 1 (i.e. {Pi{n))^'^^^ C Pi(n + 2)). Likewise, 

• The elements of (Pi(n))J-'= {k < i) begin with a stair of length exactly k — 1 and 

hence arc in Pfe(n + 1). 

From these remarks we deduce the result: 
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Figure 11: Successive copies of a Qij. 



Theorem 4 Let i be an integer such that iiiili < n + 2. Then: 

Pi{n + 2) = T{n + 2)U{Pi{n-i + l))i-^l^-UU+i u [ |J (Pfc(n + 1))^"+' J (1) 

\k>i ) 

where T(n + 2) is a set that contains one partition at most, namely: 

T{n + 2) = l {ik,k-l,...,2,l)} if 3 k integer s.t. n + 2 = 
\ otherwise 

with the initial conditions: 

■ Pi(l) = {(l,0,...)} 
=0 Vi > 1 
Pi{2) = Vi. 

Proof: Pi{n + 2) contains each of the sets in the right hand side of Equation 1. Indeed, 

if we add one grain on each of the i + 1 first columns of a sand pile with n — i + 1 
grains that begins with a stair of length i, we obtain a sand pile with n grains which 
also begins with a stair of length i. If we take a sand pile of n + 1 grains which begins 
with a stair of length at least and we add a grain on column k + 1 with k > i, we 
obtain a sand pile of n + 2 grains which also begins with a stair of length at least i. 
Finally, if n + 2 has the form for some integer k and if i is smaller than or equal 

to k (the length of the sand pile {k,k— 1, . . . , 1)) then this element of T{n + 2) begins 
with a stairs of length at least i. 

Likewise, each element of Pi{n + 2) is in one of those sets. Let s be in Pj(n + 2). 
Three cases are possible: 

• s has a step at each column, i.e. s S T(n + 2). 
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• s begins with stairs of length k with k >i and s has a plateau at fc + 1. Then, it 
is an element of {Pk+i{n + We know that such elements exist from the 
characterisation of Theorem 2. 

• s begins with a stair of length k with k > i and ,s has a cliff at k + 1. Then, s 
is an element of {Pi{n — i+ l))J-iJ-2--4i+i_ \yg know that such elements exist from 
the characterisation of Theorem 2. 

Now, let us show that the unions in Formula 1 are disjoints. The elements of 

the set {Pk{n + 1))^''+^ with k > i begin with stairs of length exactly k. So, the set 
{Pk{n + 1))^''+^ and {Pk'{n+l)y'''+^ with k,k' > i are pairwise disjoints. Moreover, the 
set {Pi{n — i + l))iiJ-2J- ■ J-iJ-i+i only contains elements that begin with stairs of length 
exactly i, so they doesn't intersect the parts (Pfc(n + l))J''=+i which begin with stairs 
of length k with k > i. Finally, if T(n + 2) is non-empty, its element clearly does not 
belong to any of the other sets. □ 

This theorem gives a better understanding of the structure of the lattices SPM{n). 

Since the unions are disjoints, the formula is even more interesting as it gives a way to 
compute the cardinality of SPM{n). We first state the following corollary, immediate 
from Theorem 4. 

Corollary 1 Let Pi^n denote |Pj(n)|, i.e. Pi^n is the number of partitions in SPM{n) 
that begin with stairs of length at least i. We have: 



Pi,n+2 = Pi,n-i+l + ^Pj,n+1 + ^i,n+2 

where 



3>i 



Si. 



n+2 



1 if3k e N s.t. n=^^^ 
otherwise 



with the initial conditions pi^i = 1, pij = for all j > 1, pj^2 = for all j > 0, and 

Po,2 = 2. 

This corollary gives a way to compute the number of elements in SPM (n) since | SPM ( 1 ) | 
1 and 

\SPM{n)\ = \SPM{n - 1)| + 

i 

Notice that this formula is nothing but the formula of the Corollary 1 specified for i = 0. 
This is not surprising, since any element of SPM{n) begins with stairs of length at least 
0. 



3 Infinite extension of SPM 

Let us now present an infinite extension of SPM. Two different possible generalisations 
are natural to extend the notions studied until here. The first one is to consider a 
column with an infinite number of grains as the initial configuration, and then study 
the evolution of the system with respect to the SPM rule. We call this model SPM{oo). 
The second one is to use the construction detailed in the first part of the paper to 
extend the order on |Jra>o '^^-^('^)- turns out that these two ideas lead to two 
isomorphic objects. This gives us an efficient way to construct SPM{oo), as shown 
below. Afterwards, we introduce the infinite tree 5PT(cx)), and we show a possible 
coding of SPM{oo) using this tree. The study of the properties of this tree gives a new 
recursive formula to compute \SPM{n)\. 
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3.1 The infinite lattice SPM (oo) 



SPM{n) is the lattice of the configurations reachable from the partition (n) by iteration 
of the SPM rule. We will now define SPM{oo) as the set of all configurations reachable 
from (oo) (this is the configuration where the first column contains infinitely many 

grains). The covering relation on SPM {oo) is defined by: s — ^ t if and only if t 
is obtained from s by application of the SPM rule on the i-th column. The order on 
SPM{oo) is the refiexive and transitive closure of this covering relation. Notice that 
any element s of SPM{oo) has the form (oo, S2, S3, . . . , Sfc)- The first partitions in 
SPMipo) are given in Figure 12 along with their covering relations (the first column, 
which always contains an infinite number of grains, is not represented on this diagram). 



1 




5 31 
6 41 22 




7 51 32 211 

8 61 42 311 




9 71 52 33 411 221 




10 81 62 43 511 321 




11 91 72 53 611 421 




12 10,1 82 63 711 44 521 331 




13 11,1 92 73 811 54 621 431 322 



Figure 12: First elements and transitions of SPM{oo). 

Notice also that the first column does not infiuence the characterisation of the ele- 
ments given in Theorem 2. We will now show that SPM{oo) is a lattice. To do so, we 
will use the notion of shot vector (see [Eri93], section 5.3). The shot vector k{s,t) from 
the sand pile s S SPM{n) to the sand pile t G SPM{n) is defined by the following: 
the i-th component ki{s,t) of k{s,t) is the number of applications of the spm rule on 
column i in order to obtain t from s. 

Wc need here an extension of this definition: the z-th component of the shot vector 
A;((oo), s) from (oo) to s G SPM{oo) is the number of applications of the spm rule on 
column i in order to obtain s from (oo). It is straightforward to see that A;((oo),s) is 
given by: 

fci((oo),s) = S2 + .S;5 + . . . 
fc2((00),s) =S3+S4 + ... 

and in general: 

fci((oo),s) = ^Sj 
j>i 

Prom [GMP98b], we have: 
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Lemma 1 Let s and t be two elements of SPM{oo). Then, 

s>t iff k{{oo),s) <k{{oo),t). 

Moreover, if m denotes the max of k{{oo), s) and k{{oo),t) then the partition u such 
that k{{oo), u) =m is in SPM{(x>) and u = inf (s, t). 

With this result, we can show that SPM{oo) is a lattice: 

Theorem 5 The set SPM{oc) is a lattice. Moreover, let s = {oo, S2, ■ ■ ■ , s^) and 
t = (00,^2, ■ ■ ■ ,ti) be two elements of SPM(oq), then, \ni{s,t) = u in SPMipo), where 

Ui = max{^^ '^i' 5^ ^i) ~ 5Z ^ ^^^'^ ^^^^ 2 < i < max{k, I) 

and sup(s,t) = inf{?i G SPM{oo),u > s,u> t}. 

Proof: From Lemma 1 and the definition of the shot vectors in SPM (00), we have the 
formula for the infimum. Since (00) is the maximal element of S'PM(oo), this set is a 
lattice. □ 

Prom the definition, it is possible to show that SPM{oo) contains an isomorphic 
copy of SPM{n) for any integer n. 

Proposition 8 Let n he a positive integer. The application: 

tt: SPM{n) — > SPM{oo) 

S = {S1,S2,- ■ ■ ,Sk) >■ S= {00,S2,...,Sk) 

is a lattice embedding, which means that it is injective and preserves the infimum and 

the supremum. 

Proof: Again, we will use the shot vector k{a,b) from a to b. Recall that ki(a,b) is 
nothing but the number of grains falling from column i in order to obtain h from a. Let 
s and t be in SPM{n). Suppose s = t, i.e. Si = ti for all i>2. We have 

si = n - ^ Si = n - ^ = ti, 

i>2 i>2 

hence s — t and tt is injective. 

It is clear that a, b in SPM{n) if and only if d — ^ b in SPM{oo), hence a >b 
in SPMin) if and only if a > 6 in SPM{oo). So, in this case, k{a,b) = k{d,b). Let 
c be the infimum of two elements a and b in SPM{n). We show that c = inf (a, 6) in 
SPM{oo). Since c = inf (a, 6) in SPM{n), we have: 

fc((n), c) = max{k{{n), a), k{{n), b)). 

Moreover, k{{n),a)) = fc((cx)),a) for all element a of SPM{n}, hence we can deduce 
fc((oo),c) = max{k{{oo),a),k{{oo),b), and c = inf(a,6) in SPM{oo), as expected. So, 
the infimum is preserved. 

Now, let us prove that the supremum is preserved. Let d = {di, . . . , d^) = sup(a, b) 
in SPM{n), and let e = (00, 62, . . . , e„i) = sup(a, b) in SPM{(x). We must show that 
d = e. To do so, we show that an element / of SPM(n) such that e = / exists. Since 
d > a and d > b in SPM{n), we have d > d and J > 6 in SPM{oo), hence d > e 
in SPM{oo), and k{{oo),d) < fc((oo),e). We first show that fci((oo),d) = fci((oo),e). 
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Without loss of generality, we can suppose that ai > bi. Notice that the partition 
(ai, ai, ai — 1, ai — 2, . . . ) is greater than a and b, hence greater than d, and so oi > di. 
Since d > a, wc have di = ai and fci((oo), d) = ki{{oo), a). Moreover, a < e < d, hence 
A;i((oo), d) = A:i((oo), e). Let us define / = (n — e2 — . . . — Cm, 62, . . . , em)- Since e < d, 
we have 62 < ^2 and since e verifies the characterisation of Theorem 2, so does /, hence 
/ G SPM{n) and e = /. Since e > a and e > 5, we have f > a and f >b, hence f > d. 
This implies that e = f > d. Therefore e = d. This gives the result. 

□ 

Let SPM{n) denotes the image by tt of SPM{n) in 5'PM(cx)). Prom Proposition 
8, SPM{n) is a sublatticc of S'PM(oo). By Theorem 4, 5'PM(n)^i is a sublattice of 
SPM{n + 1), hence, since SPM{nY^ = SPM{n), we have an increasing sequence of 
sublatticcs: 



SPM{0) < SPM{1) <■■■< SPM{n) < SPM{n +!)<•••< SPM{oo), 

where < denotes the sublatticc relation. 

Let s = (00, S2, S3, . . . , Sft) be an element of SPM{oo), then s verifies the charac- 
terisation of Theorem 2. If one takes si = S2 + 1 and n = J2i=i have that 
s' = (si, S2, . . . , Sk) is an clement of SPM{n). This implies that s = it{s') and that s is 
an element of SPM(n), therefore: 

IJ SPM{n) = SPM{oo) 

n>0 

Let us now study the disjoint union of the lattices SPM{n) for n > 0. Let us define 

S=\J SPM{n), 

n>0 

on which we extend the order relation of each SPM{n) as follows. Let s G SPM(m) 

and t € SPM{n). We define s tinSif and only if we are in one of the two following 
cases: 

• n = m and s — ^ t in SPM{n) 

• i = 0, n = m-|-l and b — a-^^ . 

In other terms, the elements of SPM{n) are linked to each other as usual whereas each 
element a of SPM{n) is linked to a^^ € SPM{n) by an edge labelled 0. Prom this 
covering relation, one can define an order on the set S as the reflexive and transitive 
closure of this covering relation. 

Theorem 6 For all integer n, SPM{n) is a sublattice of S. 

Proof: The fact that SPM{n) is present in S is immediate from the definition. What 
we have to show is that the lattice structure of SPM{n) is preserved in S. Let s and t 
be two elements of SPM(n). We have to show that the infimum and the supremum of 
s and t in S are in SPM(n). Let u = inf(s,i) in SPM{n) and u' = inf(s,t) in S. We 
have that s > u' > u in S, thus 

i>l i>l i>l 
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and 

i>l 

Therefore u' is an element of SPM{n), and we have u' = u. The same method can be 
applied for the supremurn. q 

This result is illustrated in Figure 13 (right). The surprising result is that these two 
ways to extend the sand pile model to infinity, i.e. the first one by adding new elements 

to some SPM{n) to extend it into SPM{n + 1) and infinitely iterating the process to 
obtain SPM{oo), and the second one by linking together all the SPM(n) for all n to 
obtain S, lead to the same object. 

Theorem 7 The application x defined by: 

X:S=\_\ SPM{n) — > SPM{oo) 

n>0 

S = (Sl, S2, • • • , Sfe) 1-^ X(s) = (OO, Si, S2, . . . , Sk) 

is a lattice isomorphism, which means that it is one to one and preserves the supremum 

i i I 1 

and the infimum. Moreover, s — > t in |Jn>o SPM(n) if and only if x(s) — »■ x(t) in 

SPM{oo). 

Proof: The application x is obviously injective. Let us show that x is surjective. 
Let s = (cx),si,... ,si) be an element of SPM{oo). Define n = J2i>i^i ^-i^^ — 
(si, . . . , si). Since s is in SPM{oo), s verifies the conditions of Theorem 2, and so does 
s'. Therefore s' is an element of SPM(n), and since s = x(s')) the application x is 
surjective. 

It is clear that for all s,t G S, one has s — ^ t if and only if x(s) - — > x{t)- Obviously, 
X is an order isomorphism. Since SPM{oo) is a lattice, we can conclude that X is a 
lattice isomorphism. q 

This result is illustrated in Figure 13 (left). In the following, we simplify the nota- 
tions by representing the elements of S'PM(cx)) without their first column. Our aim is 
now to construct large parts of SPM{oo). A first solution is to construct SPM{n) for 
large values of n. However, this does not lead to filters (a filter of a lattice is a subset 
of this lattice closed for the supremum) of SPM{oo). We will now define special filters 
of SPM((x) and explain how we can construct them efficiently. For a given n, let us 
denote by SPM{< n) the set 

IJ SPM{i). 

0<i<n 

For example, SPM{< 7) is shown in Figure 13. It is easy to see that SPM{< n) is a 
filter of SPM{oo) for all n. The infinite lattice SPM{oo) can be regarded as a limit of 
this sequence of poscts. The results presented in this section give us an efficient method 
to construct SPM{< n) for all n (see Algorithm 2 and Theorem 8). Moreover, we show 

another property of SPM{< n). 

Proposition 9 The poset SPM{< n) is a sublattice of SPM{<x>) for all n. 

Proof: To show the claim, it suffices to consider s G SPM{k) and t 6 SPM{1), with 
k < I < n, and show that inf(s,i) and sup(s,t) (which are in SPM{oo) since SPM (oo) 
is a lattice) are also in SPM{< n). Let u = sup(s,i) in SPM{oo). Since u > s in 
SPM{oo), there exists an integer m < k such that u £ SPM{m), therefore u is in 
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SPM{< n). Let now u' = mf{s,t) in SPM{oo). Let s' be the partition obtained by 
addition of I — k grains on the first column of s. Then, s' S SPM{1) and inf (s', t) in 
SPM{oo) is also in SPM{1). So, u = mi{s,t) is also greater than or equal to inf(s',f) 
in S'PM(cx)), hence u e S'PM(< /). □ 



Algorithm 2 CONSTRUCTION OF SPM{< n) 
Input: an integer n 
Output: SPM{< n) 
begin 

- Init Result with S'PM(O); 

- for 1 < i < n do 

- Extract SPM{i) from Result: 

Start from (0); 

Depthfirst search Result to get the connex part containing SPM{i); 

- Link SPM{i) to Result: 
for each s e 5'PM(i - 1) do 

|_ Link s G SPM{i - 1) to s^^ m SPM{i) with an edge -U; 

- return Result 
end 




Figure 13: The two ways to identify SPM{n) in SPM{oo) for all n. 



Theorem 8 Algorithm 2 computes SPM{< n) in time linear in the number of elements 
and edges in SPM{< n). 
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Proof: Let l^l denote the number of elements and edges in S for any lattice S. For 
all i between 1 and n, the lattice SPM{i) is a connected component of Result. It is 
well known that the computation of a connected component has a cost proportional 
to the cardinality of this same component, so its extraction can be processed in linear 
time. The addition of the edges that link SPM{i) to SPM{i — 1) can be done in 
0{\SPM{i — 1)1). Therefore, each iteration of the for loop is executed in 0{\SPM{i)\)^ 
hence the execution of the whole loop is linear with respect to the total number of 
elements and edges of SPM{< n). This is the asymptotic cost of the whole algorithm. 

□ 

3.2 The infinite tree SPT{oo) 

As shown in our construction oi SPM(n+l) from SPM{n), each element s of SPM(n+ 
1) is obtained from an element s' G SPM{n) by addition of one grain: s = s'^* with 
i an integer between 1 and e{s') + 1. Thus, we can define an infinite tree SPT{oo) 
(for Sand Pile Tree) whose nodes are the elements of |Jn>o SPM{n) and in which the 
fatherhood relation is defined by: 

t is the i-th son of s if and only if t = s^* for some i, 1 <i < e(s) + 1. 

The edge s — > s^' is labelled with i. The root of this tree is (0). The eight first levels 

of SPT{oc) arc shown in Figure 14 (we call the set of elements of depth n the "level n" 
of the tree). Each node s of SPT{oo) has e(s) + 1 sons linked to s with edges labelled 

1 2 e(s) + l 

. . . , >. 

Notice that, although the notation is the same as the one used for the spm transitions 
in the lattice, an edge s t in the tree means that t is obtained from s by addition 
of one grain on its i-th column {t = s-'-*), and not that t is obtained from s by having 
the top grain of the i-th column fall onto the {i + l)-th. Therefore, if s € SPM{n) then 
t £ SPM{n + 1). So, the structure of the lattices is not directly visible in SPT{oo). 
One goal in the following part of the section will be to explore the possibility of the 
construction of the lattices from the tree. The first results will be to state, as we could 
have guess from Theorem 7, that there are two ways to find SPM{n) in SPT{oo). 
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Figure 14: The first levels of SPT{oo). 

Proposition 10 The level n of SPT (oo) contains exactly the elements of SPM{n). 

Proof: straightforward from the construction of SPM{n + 1) from SPM{n) given 
above. □ 

21 



Proposition 11 For all integer n, the set SPM{n) is a subtree of SPT {oo) that con- 
tains its root. 

Proof: The proposition is obviously true for n = 1. Let us suppose it is true for n, and 
let us verify it is true for n + 1 . By construction, the elements of SPM (n + 1 ) \ SPM (n) 
are sons of elements of SPM{n). The result follows. q 

We will now show that SPT{oo) can be described recursively, which allows us to give 
a new recursive formula for \SPM{n)\. Let us first consider one element s of SPT{oo) 
and let k = e(s). By definition of SPT{oo), s has exactly k + 1 sons. Notice that the k 
first sons of s all verify the same following property: for all integer i between 1 and k, 
the son s^' begins with a stair of length i — 2, has a plateau at z — 1, and a cliff at i. 
Prom this remark, we introduce certain types of subtrees of SPT{oo). 

Definition 1 We call N). subtree, with k >1, any subtree T of SPT {oo) which is rooted 
at an element s that begins with a stair of length k — 2, has a plateau at k — 1 and a 
cliff at k, and contains all the descendants of s. 

Definition 2 We call Xk subtree, with k > 1, any subtree T of SPT (oo) which is 
rooted at a node that has at least k sons and such that the i-th son is the root of a Ni 
subtree for all i between 1 and k. The structure of Xk subtrees is shown in Figure 15. 
Moreover, we define Xq = %. 




Figure 15: Definition of the Xk subtrees. 

Notice that SPT{oo) contains some Xk subtrees. Indeed, if one takes an element s 
which begins with a stair of length k, then its i-th son is the root of a subtree for all 
i, 1 < i < k. So s is the root of a Xk subtree. Notice also that s might have other sons 
outside the Xk subtree. With this remark, we characterise in the following proposition 
the structure of the Nk subtrees, with fc > 1. See Figure 18 and 19. 

Proposition 12 We have the following statements : 

1. A Ni subtree is an infinite chain whose edges are all labelled 1. 

2. A Nk subtree, with k > 2, is com,posed by a chain of k nodes whose edges are 
labelled k — 1, k — 2, . . . , 1 and whose i-th node is the root of a Xk-i-i subtree 
for all i between 1 and k — 2. Moreover, the k-th node is root of a Xk subtree. 

Proof: 

1. By definition, the root of a A^i subtree is an element s which begins with a cliff, 
e(s) + 1 = 1 (the conditions on the initial stair and the plateau make no sense in 

this case). Therefore, the only son of s is s — s^^, which also begins with a cliff. 
Therefore, s^^ is the root of a A^i subtree and we can conclude by induction. 
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2. Let k be greater than 1. Let us consider a partition s such that: 



• e(s) = fc - 2, 

• s has a plateau at fc — 1, 

• and s has a chfF at k. 

The node s has A; — 1 sons: s^^ , s-^^, . . . , s^''-^. Since e(s) = A; — 2, and from 
the remarks above, the node s is the root of a Xk-2 subtree that contains all the 
elements reachable from its fc — 2 first sons. Let t be the (fc — l)-th son of s (t is 
outside the Xk-2 subtree rooted at s). The subtree rooted at s is then the union 
of a Xk-2 subtree and a subtree with root t (see Figure 16). Look now at the 




Figure 16: First step of the structure of Nk- 

subtree with root t. Wo have that e{t) = fc — 3, hence t is the root of a Xk-3 
subtree. Let u be the (fc — 2)-th son of t. We then obtain Figure 17. When this 




Figure 17: Second step of the structure of Nk- 

process is iterated, we obtain x = si'=-iifc-2 -42_ 'pj-^jg element x begins with a 
plateau of length 1 followed by a stair of length fc — 2 and a cliff. Therefore x has 
only one son x — > y. This element y begins with a stair of length fc — 1 followed 
by a cliff. As noticed above, y is the root of a Xk subtree and this subtree contains 
all the elements reachable from y. Then we obtain the announced structure of Nk 
subtrees, see Figure 19. q 

Figure 18: Structure of the Ni subtrees. 

Using now the fact that a Xk subtree is defined in terms of Nk subtrees, we can 
describe the structure of a Nk subtree only in terms of other N^ subtrees with i < fc as 
shown Figure 20. 

Notice that we can deduce directly by induction from the structure of the Nk subtrees 
shown in Figure 20 that all the edges in a A''^ subtree are labelled with integer smaller 
than or equal to fc. The recursive structures we have defined, and the propositions given 
above allow a compact representation of the tree SPT{oo) as a chain: 

Theorem 9 The tree SPT{oo) can he represented by the infinite chain shown in Figure 
21. The nodes of this chain are the fixed points of SPM{n) for n > 0. 
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Figure 19: Structure of Nk subtrees. 




Figure 20: Self referencing structure of Nk- 



The chain is defined as follows: let k a positive integer and let Pk = {k,k — l,k — 
2, . . . , 2, 1) and Pk+i = {k + l,k,k — l,k — 2, . . . ,2, 1); the subchain between Pk and 
Pk+i contains k + 1 nodes: 

P ^~^] plk+l , plk+llk k-1 2 ^ plk+l---l2 1 ^ p ^ 

where each node p^'^+i '''-' with i between 3 and k + 1 is the root of a Xi-2 subtree, and 
Pk is the root of a Xk subtree. 

Proof: Let us consider the rightmost chain in SPT{oo). This chain is composed by 
the fixed points of SPM(n) for n > 0. Let fc be a positive integer. Let us consider 
the subchain of this chain that begins with Pk = (fc, fc — 1, . . . ,1) and terminates with 
Pfe+i = (fc+l,fc,... ,1): 

rj fe+\ r>lk+i k r,lk+ilk k-\ 2 r>lk+i---l2 1 p 

n — > I'k — * ^k — — ^ ^k — * n+i- 

The node s = p^''+'^---^' with i >2, begins with stairs of length i — 2 followed by a 
plateau at i — 1, hence s is the root of a Xi^2 subtree and its last son is obtained by 
plk+i-li izi plk+i-li-i _ rpj^jg jg ^jjg j^g^^ jjQjjg ^jjg ^.j^^jjj Therefore SPT{oo) can 
be described as indicated. □ 

As seen above, the level n of SPT{oo) contains exactly SPM{n). Therefore, it 
suffices to count the number of paths of length n from the root of SPT{oo) to obtain 
\SPM{n)\. The recursive structure of the tree, detailed above, gives us a way to achieve 
this. 
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11 21 211 221 321 3211 3221 






Figure 21: The tree SPT{oo) represented as an infinite chain. 

Theorem 10 Let c{l, k) denote the number of paths in a subtree originating from 
the root and of length I, then we have: 



c{l, k) = < 



ifl<Oandk<0 

1 if I > and k = 1 
k ifl = l and k > 

c{l — k,k) + '^'^~^ c{l — i + l,k — i) + e{l,k) otherwise 



where e{l, k) = if k > I and e{l, k) = 1 otherwise. 

Proof: The proof follows from the recursive structure of Xk detailed above. There is 
no path of length or less, and the Xk are empty for fc < 0, hence the first case. A 
Xi subtree is a simple chain, hence there is exactly 1 path of any length, hence the 
second case. The third case is immediately deduced from the fact that the root of a Xk 
subtree has exactly k sons. Finally, the recursive formula in the fourth case comes from 
the fact that the structure of Xk subtrees shown in Figure 15 allows us to consider a 
Xk subtree as a node s where s is the root of a Xk-i subtree and has one more son 
which is the root of a Nk subtree. Then, from the structure of Nk subtrees in terms 
of Xk ones shown in Figure 19, we deduce a description of Xk subtrees in terms of Xi 
subtrees with <i < k, from which the formula is straightforward. q 

Corollary 2 The cardinality of SPM[n) is given by: 

|5PM(n)| = l + ^^c(n-^^^-j + l,i-i + l) 
where k is the integer such that ^s!Js+ll <n< (^±iMM. 

Proof: This formula is deduced from the chain structure of the tree, shown in Figure 
21. The quantity 1 corresponds to the path of length n that follows the chain without 
entering in a Xi subtree. The double sum corresponds to the repartition of the Xi 
subtrees along the chain. q 

We will now show how information on SPM{oo) can be deduced from SPT{oo). 
The lattice structure SPM{oo) and the infinite tree SPT{oo) are defined over the same 
underlying set: Un>o '^-'^-^(^)- Therefore we can easily give a bijection from one to 
the other. We now show how the ordered structure of SPM {oo) can be deduced from 
SPT{oo). 

Proposition 13 Every element of SPT{oo) has an outgoing edge with label 1 in SPM{oo) 

Moreover, for i > 1, a partition s of SPT{oo) has an outgoing edge — > in SPM{oo) 
if and only if s belongs to a Ni subtree of SPT{oo). It should be noticed that s may 
belong to several Ni subtrees for dictinct values of i. In this case, s is the origin of an 

edge in SPM{oo) for each such i. In other words, a is in a Ni subtree of SPT{oo) 
if and only if there is an outgoing edge from a labelled with i + 1 in SPM{oo). 
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Proof: By definition of S'PAf (oo), every partition in SPAI{oo) obviously has a succes- 
sor by — Suppose now i > 1 and suppose that a is in a Ni subtree of SPT{oo), and 
denote by s the root of this Ni subtree. By definition of the Ni subtrees, we have that 
Si — Si+i > 2, hence — a^+i > 2. Therefore, when adding the first infinite column 
to obtain the corresponding element of SPM{oo), we obtain that there is an ouot going 
edge from a labelled i + 1 in SPM{oo). 

Let us now consider a partition a having an outgoing edge labelled with ? + 1 in 
SPM(oo). From the structure of the tree shown in Figure 21, we know that a is in 
a Xk subtree of SPT{oo). Indeed, the elements of the chain arc fixed points of SPM, 
hence they have no outgoing edge in SPM{oo) except the ones labelled 1, as said above. 

Let s be the root of such a subtree of SPT{oo), that is, s is the root of a 
subtree that contains a. Consider a path from s to a in this tree: 

ii i2 is 'i-i ii 

s = So — >■ Si — »■ S2 — >• . . . — > Si_i — > Si = a. 

We have that Si — Si+i < 1 by definition and ai — Ui+i > 2 since there is a transition 
from a labelled i + 1 in SPM{oo). Therefore, there exists an integer j such that ij — i. 
In fact, we have an even stronger condition on the path: it must verify |« + 1| < \i\ 
where \x\ denotes the number of edges labelled with x on the path. But it is easy to 
see from the structure of the Ni subtrees shown in Figure 20 that the only case where 

this happends is when Sj sj+i such that sj+i is the root of a Ni subtree, hence a 
is in a Ni subtree, as announced. q 

It follows from this proposition that we can find all the (immediate) successors of a 
partition s in SPM{(X>): it is sufficient to go from the root of SPT (oo) to s G SPT (oc) 
and so determinate the integers i such that s is in a Ni subtree. 



4 Conclusion and Perspectives 

Through the study of the construction of SPM(n+l) from SPM{n), we obtained much 
information about this set. First, it is strongly self-similar and can be constructed 
using this property. Moreover, this construction procedure gives a formula for the 
cardinal of SPM{n), where no formula was known before. In a second part, we gave a 
natural way to extend SPM(n) to infinity, and again self-similarity of this infinite lattice 
appeared. Finally, we gave a tree structure to the sets SPM{n) and SPM{oo), which 
allows efficient enumeration of SPM{n), as well as another formula for the cardinal of 
SPM{n). 

The duplication process that appears during the construction of the lattices SPM(n) 
may be much more general, and could be extended to other kinds of lattices, maybe 
leading to the definition of a special class of lattices, which contains the lattices SPM(n). 
Moreover, the ideas dcvelopped in this paper could be applied to others dynamical 
systems, such as the Brylawski dynamical system [Bry73], Chip Firing Games, or tilings 
with flips, with some benefit. 
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